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METHODS FOR MOTION ESTIMATION WITH 
ADAPTIVE MOTION ACCURACY 



BACKGROUND OF THE INVENTION 
5 The present invention relates generally to a method of compressing or 

coding digital video with bits and, specifically, to an effective method for estimating and 
encoding motion vectors in motion-compensated video coding. 

In classical motion estimation the current frame to be encoded is 
decomposed into image blocks of the same size, typically blocks of 16x16 pixels, called 

1 0 "macroblocks." For each current macroblock, the encoder searches for the block in a 
previously encoded frame (the "reference frame") that best matches the current 
macroblock. The coordinate shift between a current macroblock and its best match in the 
reference frame is represented by a two-dimensional vector (the "motion vector") of the 
macroblock. Each component of the motion vector is measured in pixel units. 

1 5 For example, if the best match for a current macroblock happens to be at 

the same location, as is the typical case in stationary background, the motion vector for 
the current macroblock is (0,0). If the best match is found two pixels to the right and 
three pixels up from the coordinates of the current macroblock, the motion vector is (2,3). 
Such motion vectors are said to have integer pixel (or "integer-pel" or "full-pel") 

20 accuracy, since their horizontal X and vertical Y components are integer pixel values. In 
FIG. 1, the vector V t =(1,1) represents the full-pel motion vector for a given current 
macroblock. 

Moving objects in a video scene do not move in integer pixel increments 
from frame to frame. True motion can take any real value along the X and Y directions. 
25 Consequently, a better match for a current macroblock can often be found by 

interpolating the previous frame by a factor NxN and then searching for the best match in 
the interpolated frame. The motion vectors can then take values in increments of 1/N 
pixel along X and Y and are said to have 1/N pixel (or "1/N-pel") accuracy. 
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In "Response to Call for Proposals for H.26L," ITU-Telecommunications 
Standardization Sector, Q.15/SG16, doc. Q15-F-11, Seoul, Nov. 98, and "Enhancement 
of the Telenor proposal for H.26L," ITU-Telecommunications Standardization Sector, 
Q.15/SG16, doc. Q15-G-25, Monterey, Feb. 99, Gisle Bjontegaard proposed using 1/3- 
5 pel accurate motion vectors and cubic-like interpolation for the H26L video coding 
standard (the "Telenor encoder"). To do this, the Telenor encoder interpolates or "up- 
samples" the reference frame by 3x3 using a cubic-like interpolation filter. This 
interpolated version requires nine times more memory than the reference frame. At a 
given macroblock, the Telenor encoder estimates the best motion vector in two steps: the 

10 encoder first searches for the best integer-pel vector and then the Telenor encoder 

searches for the best 1/3-pixel accurate vector V 1/3 near Using FIG. 1 as an example, 
a total of eight blocks (of 16x16 pixels) in the 3x3 interpolated reference frame are 
checked to find the best match which, as shown is the block associated to the motion 
vector V 1/3 = (VX, VY) = (1+1/3,1). The Telenor encoder has several problems. First, it 

15 uses a sub-optimal fast-search strategy and a complex cubic filter (at all stages) to 

compute the 1/3-pel accurate motion vectors. As a result, the computed motion vectors 
are not optimal and the memory and computation requirements are very expensive. 
Further, the Telenor encoder uses an accuracy of the effective rate-distortion criteria that 
is fixed at 1/3-pixel and, therefore, does not adapt to select better motion accuracies. 

20 Similarly, the Telenor encoder variable-length code ("VLC") table has an accuracy fixed 
at 1/3-pixel and, therefore, is not adapted and interpreted differently for different 
accuracies. 

Most known video compression methods estimate and encode motion 
vectors with 1/2 -pixel accuracy, because early studies suggested that higher or adaptive 
25 motion accuracies would increase computational complexity without providing additional 
compression gains. These early studies, however, did not estimate the motion vectors 
using optimized rate-distortion criteria, did not exploit the convexity properties of such 
criteria to reduce computational complexity, and did not use effective strategies to encode 
the motion vectors and their accuracies. 
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One such early study was Bernd Girod's "Motion-Compensating 
Prediction with Fractional-Pel Accuracy IEEE Transactions on Communications, 
Vol. 41, No. 4, pp. 604-612, April 1993 (the "Girod work"). The Girod work is the first 
fundamental analysis on the benefits of using sub-pixel motion accuracy for video 
5 coding. Girod used a simple, hierarchical strategy to search for the best motion vector in 
sub-pixel space. He also used simple mean absolute difference ("MAD") criteria to select 
the best motion vector for a given accuracy. The best accuracy was selected using a 
formula that is not useful in practice since it is based on idealized assumptions, is very 
complex, and restricts all motion vectors to have the same accuracy within a frame. 

1 0 Finally, Girod focused only on prediction error energy and did not address how to use 
bits to encode the motion vectors. 

Another early study was Smita Gupta's and Allen Gersho's "On Fractional 
Pixel Motion Estimation," Proc. SPIE VCIP, Vol. 2094, pp. 408-419, Cambridge, 
November 1993 (the "Gupta work"). The Gupta work presented a method for computing, 

15 selecting, and encoding motion vectors with sub-pixel accuracy for video compression. 
The Gupta work disclosed a formula based on mean squared error ("MSE") and bilinear 
interpolation, used this formula to find an ideal motion vector, and then quantized such 
vector to the desired motion accuracy. The best motion vector for a given accuracy was 
found using the sub-optimal MSE criteria and the best accuracy was selected using the 

20 largest decrease in difference energy per distortion bit, which is a greedy (sub-optimal) 
criteria. A given motion vector was coded by first encoding that vector with 1/2-pel 
accuracy and then encoding the higher accuracy with refinement bits. Course-to-fine 
coding tends to require significant bit overhead. 

In "On the Optimal Motion Vector Accuracy for Block-Based Motion- 

25 Compensated Video Coders," Proc. IST/SPIE Digital Video Compression: Algorithms 
and Technologies, pp. 302-314, San Jose, February 1996 (the "Ribas work"), Jordi Ribas- 
Corbera and David L. Neuhoff, modeled the effect of motion accuracy on bit rate and 
proposed several methods to estimate the optimal accuracies that minimize bit rate. The 
Ribas work set forth a full-search approach for computing motion vectors for a given 



accuracy and considered only bilinear interpolation. The best motion vector was found 
by minimizing MSE and the best accuracy was selected using some formulas derived 
from a rate-distortion optimization. The motion vectors and accuracies were encoded 
with frame-adaptive entropy coders, which are complex to implement in real-time 
5 applications. 

In "Proposal for a new core experiment on prediction enhancement at 
higher bitrates," ISO/IEC JTC1/SC29/WG1 1 Coding of Moving Pictures and Audio, 
MPEG 97/1827, Sevilla, Feb. 1997 and "Performance Evaluation of a Reduced 
Complexity Implementation for Quarter Pel Motion Compensation," ISO/IEC 

10 JTC1/SC29/WG1 1 Coding of Moving Pictures and Audio, MPEG 97/3 146, San Jose, 
Jan. 1998, Ulrich Benzler proposed using 1/4-pel accurate motion vectors for the video 
sequence and more advanced interpolation filters for the MPEG4 video coding standard. 
Benzler, however, used the Girod's fast-search technique to find the 1/4-pel motion 
vectors. Benzler did consider different interpolation filters, but proposed a complex filter 

15 at the first stage and a simpler filter at the second stage and interpolated one macroblock 
at a time. This approach does not require much cache memory, but it is computationally 
expensive because of its complexity and because all motion vectors are computed with 
1/4-pel accuracy for all the possible modes in a macroblock (e.g., 16x16, four-8x8, 
sixteen-4x4, etc.) and then the best mode is determined. Benzler used the MAD criteria 

20 to find the best motion vector which was fixed to 1/4-pel accuracy for the whole 

sequence, and hence he did not address how to select the best motion accuracy. Finally, 
Benzler encoded the motion vectors with a variable-length code ("VLC") table that could 
be used for encoding 1/2 and 1/4 pixel accurate vectors. 

The references discussed above do not estimate the motion vectors using 

25 optimized rate-distortion criteria and do not exploit the convexity properties of such 
criteria to reduce computational complexity. Further, these references do not use 
effective strategies to encode motion vectors and their accuracies. 
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BRIEF SUMMARY OF THE INVENTION 

One preferred embodiment of the present invention addresses the 
problems of the prior art by computing motion vectors of high pixel accuracy (also 
denoted as "fractional" or "sub-pixel 11 accuracy) with a minor increase in computation. 
5 Experiments have demonstrated that, by using the search strategy of the 

present invention, a video encoder can achieve significant compression gains (e.g., up to 
thirty percent in bit rate savings over the classical choices of motion accuracy) using 
similar levels of computation. Since the motion accuracies are adaptively computed and 
selected, the present invention may be described as adaptive motion accuracy ("AMA"). 

1 0 One preferred embodiment of the present invention uses fast-search 

strategies in sub-pixel space that smartly searches for the best motion vectors. This 
technique estimates motion vectors in motion-compensated video coding by finding a 
best motion vector for a macroblock. The first step is searching a first set of motion 
vector candidates in a grid of sub-pixel resolution of a predetermined square radius 

15 centered on V\ to find a best motion vector V 2 . Next, a second set of motion vector 

candidates in a grid of sub-pixel resolution of a predetermined square radius centered on 
V 2 is searched to find a best motion vector V 3 . Then, a third set of motion vector 
candidates in a grid of sub-pixel resolution of a predetermined square radius centered on 
V 3 is searched to find the best motion vector of the macroblock. 

20 In an alternate preferred embodiment the present invention a technique for 

estimating high-accurate motion vectors may use different interpolation filters at different 
stages in order to reduce computational complexity. 

Another alternate preferred embodiment of the present invention selects 
the best vectors and accuracies in a rate-distortion ("RD") sense. This embodiment uses 

25 rate-distortion criteria that adapts according to the different motion accuracies to 
determine both the best motion vectors and the best motion accuracies. 

Still further, another alternate preferred embodiment of the present 
invention encodes the motion vector and accuracies with an effective VLC approach. 
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This technique uses a VLC table that is interpreted differently at different coding units, 
according to the associated motion vector accuracy. 

The foregoing and other objectives, features, and advantages of the 
invention will be more readily understood upon consideration of the following detailed 
5 description of the invention, taken in conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING 

FIG. 1 is a diagram of an exemplary full-pel and 1/3-pel locations in 

velocity space. 

10 FIG. 2 is a flowchart illustrating a prior art method for estimating the best 

motion vector. 

FIG. 3 is a diagram of an exemplary location of motion vector candidates 
for full-search in sub-pixel velocity space. 

FIG. 4 is a flowchart illustrating a full-search preferred embodiment of the 
1 5 method for estimating the best motion vector of the present invention. 

FIG. 5 is a diagram of an exemplary location of motion vector candidates 
for fast-search in sub-pixel velocity space. 

FIG. 6 is a flowchart illustrating a fast-search preferred embodiment of the 
method for estimating the best motion vector of the present invention. 
20 FIG. 7 is a detail flowchart illustrating an alternate preferred embodiment 

of step 114 of FIG. 6. 

FIG. 8 is a graphical representation of experimental performance results of 
the Telenor encoder with and without AMA in the "Container" video sequence, with 
QCIF resolution, and at the frame rate of 10 frames per second. 
25 FIG. 9 is a graphical representation of experimental performance results of 

the Telenor encoder with and without AMA in the "News" video sequence, with QCIF 
resolution, and at the frame rate of 10 frames per second. 
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FIG. 10 is a graphical representation of experimental performance results 
of the Telenor encoder with and without AM A in the "Mobile" video sequence, with 
QCIF resolution, and at the frame rate of 10 frames per second. 

FIG. 1 1 is a graphical representation of experimental performance results 
5 of the Telenor encoder with and without AMA in the "Garden" video sequence, with SIF 
resolution, and at the frame rate of 15 frames per second. 

FIG. 12 is a graphical representation of experimental performance results 
of the Telenor encoder with and without AMA in the "Garden" video sequence, with 
QCIF resolution, and at the frame rate of 15 frames per second. 
10 FIG. 13 is a graphical representation of experimental performance results 

of the Telenor encoder with and without AMA in the "Tempete" video sequence, with 
SIF resolution, and at the frame rate of 15 frames per second. 

FIG. 14 is a graphical representation of experimental performance results 
of the Telenor encoder with and without AMA in the "Tempete" video sequence, with 
15 QCIF resolution, and at the frame rate of 15 frames per second. 

FIG. 15 is a graphical representation of experimental performance results 
of the Telenor encoder with and without AMA in the "Paris shaked" video sequence, with 
QCIF resolution, and at the frame rate of 10 frames per second. 

FIG. 16 is a graphical representation of experimental performance results 
20 of fast-search ("Telenor FS AMA+c") and full-search ("Telenor AMA+c") strategies in 
the "Mobile" video sequence, with QCIF resolution, and at the frame rate of 10 frames 
per second. 

FIG. 17 is a graphical representation of experimental performance results 
of fast-search ("Telenor FSAMA+c") and full-search ("Telenor AMA+c") strategies in 
25 the "Container" video sequence, with QCIF resolution, and at the frame rate of 10 frames 
per second. 

FIG. 18 is a graphical representation of experimental performance results 
of tests using only one reference frame for motion compensation as compared to tests 
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using multiple reference frames for motion compensation the in the "Mobile" video 
sequence, with QCIF resolution, and at the frame rate of 10 frames per second. 

DETAILED DESCRIPTION OF THE INVENTION 
5 The methods of the present invention are described herein in terms of the 

motion accuracy being modified at each image block. These methods, however, may be 
applied when the accuracy is fixed for the whole sequence or modified on a frame-by-frame 
basis. The present invention is also described as using Telenor's video encoders (and 
particularly the Telenor encoder) as described in the Background of the Invention. Although 

10 described in terms of Telenor's video encoders, the techniques described herein are 
applicable to any other motion-compensated video coder. 

Most video coders use motion vectors with half pixel (or ,r l/2-pel") accuracy 
and bilinear interpolation. The first version of Telenor's encoder also used 1/2-pel motion 
vectors and bilinear interpolation. The latest version of Telenor's encoder, however, 

15 incorporated 1/3-pel vectors and cubic-like interpolation because of the additional 
compression gains. Specifically, at a given macroblock, Telenor's estimates the best motion 
vector in two steps shown in FIG. 2. First, the Telenor encoder searches for the best integer- 
pel vector Vj (FIG. 1) 100. Second, the Telenor encoder searches for the best 1/3-pixel 
accurate vector V 1/3 (FIG. 1) near V 1 102. This second step is shown graphically in FIG. 1 

20 where a total of eight blocks (each having an array of 16x16 pixels) in the 3x3 interpolated 
reference frame are checked to find the best match. The motion vectors for these eight 
blocks are represented by the eight solid dots in the grid centered on Y v In FIG. 1 the best 
match is the block associated to the motion vector V 1/3 = (V x , V Y ) = (1+1/3, 1). 

The technology of the present invention allows the encoder to choose 

25 between any set of motion accuracies (for example, 1/2, 1/3, and 1/6-pel accurate motion 
vectors) using either a full search strategy or a fast search strategy. 
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Full-Search AMA Search Strategy 

As shown in FIGS. 3 and 4, in the full-search adaptive motion accuracy 
("AMA") search strategy the encoder searches all the motion vector candidates in a grid 
of 1/6-pixel resolution and a "square radius" (defined herein as a square block defined by 
5 a number of pixels up, a number of pixels down, and a number of pixels to both sides) of 
five pixels as shown in FIG. 3, FIG. 4 shows that the first step of the full-search AMA is 
to search for the best integer-pel vector V x (FIG. 1) 104. In the second step of the full- 
search AMA, the encoder searches for the best 1/6-pixel accurate vector V 1/6 (FIG. 3) near 
Y l 106. In other words, the full-search AMA modifies the second step of the Telenor's 

10 process so that the encoder also searches for motion vector candidates in other sub-pixel 
locations in the velocity space. The objective is to find the best motion vector in the grid, 
i.e., the vector that points to the block (in the interpolated reference frame) that best 
matches the current macroblock. Although the full-search strategy is computationally 
complex since it searches 120 sub-pixel candidates, it shows the full potential of this 

1 5 preferred method of the present invention. 

A critical issue in the motion vector search is the choice of a measure or 
criterion for establishing which block is the best match for the given macroblock. In 
practice, most methods use either the mean squared error ("MSE") or mean absolute 
difference ("MAD") criteria. The MSE between two blocks consists of subtracting the 

20 pixel values of the two blocks, squaring the pixel differences, and then taking the 

average. The MAD difference between two blocks is a similar distortion measure, except 
that the absolute value of the pixel differences is computed instead of the squares. If two 
image blocks are similar to each other, the MSE and MAD values will be small. If, 
however, the image blocks are dissimilar, these values will be large. Hence, typical video 

25 coders find the best match for a macroblock by selecting the motion vector that produces 
either the smallest MSE or the smallest MAD. In other words, the block associated to the 
best motion vector is the one closest to the given macroblock in an MSE or MAD sense. 

Unfortunately, the MSE and MAD distortion measures do not take into 
account the cost in bits of actually encoding the vector. For example, a given motion 



9 



vector may minimize the MSE, but it may be very costly to encode with bits, so it may 
not be the best choice from an coding standpoint. 

To deal with this, advanced encoders such as those described by Telenor 
use rate-distortion ("RD") criteria of the type "distortion + L*Bits" to select the best 
5 motion vector. The value of "distortion" is typically the MSE or MAD, "L" is a constant 
that depends on the compression level (i.e., the quantization step size), and "Bits" is the 
number of bits required to code the motion vector. In general, any RD criteria of this 
type would work with the present invention. However, in the present invention "Bits" 
include the bits needed for encoding the vector and those for encoding the accuracy of the 
10 vector. In fact, some candidates can have several "Bits" values, because they can have 

several accuracy modes. For example, the candidate at location (1/2, -1/2) can be thought 
of having 1/2 or 1/6 pixel accuracy. 

Fast-Search AMA Search Strategy 

15 As shown in FIGS. 5 and 6, in the fast-search adaptive motion accuracy 

("AMA") search strategy the encoder checks only a small set of the motion vector 
candidates. In the first step of the fast-search AMA, the encoder checks the eight motion 
vector candidates in a grid of 1/2-pixel resolution of square radius 1, which is centered on 
V l 108. V 2 is then set to denote the candidate that has the smallest RD cost (i.e., the best of 

20 the eight previous vectors and 110. Next, the encoder checks the eight motion vector 
locations in a grid of 1/6-pixel resolution of square radius 1 that is now centered on V 2 1 12, 
If V 2 has the smallest RD cost 1 14, the encoder stops its search and selects V 2 as the motion 
vector for the block. Otherwise, V 3 is set to denote the best motion vector of the eight 116. 
The encoder then searches for a new motion vector candidates in the grid of 1/6-pixel 

25 resolution of square radius 1 that is centered on V 3 118. It should be noted that some of the 
candidates in this grid have already been tested and can be skipped. The candidate with the 
smallest RD cost in this last step is selected as the motion vector for the block 120. 

Experimental data has shown that, on average, this simple fast search strategy 
typically checks the RD cost of about eighteen locations in sub-pixel space (ten more than 

10 



Telenor's search strategy), and hence the overall computational complexity is only 
moderately increased. 

The experimental data discussed below in connection with FIGS, 8-18 show 
that there is practically no loss in compression performance from using this fast-search 
5 version of AMA. This is because the fast-search AMA search strategy exploits the convexity 
of the "distortion + L*Bits" curve (c.£ , "distortion" is known to be convex), by creating a 
path that smartly follows the RD cost from higher to lower levels. 

Alternate embodiments of the invention replace one or more of the steps 108- 
120. These embodiments have also been effective and have further reduced the number of 

10 motion vector candidates to check in the sub-pixel velocity space. 

FIG. 7, for example, checks candidates of 1/3-pel accuracy. In this 
embodiment step 1 12 is replaced by one of three possible scenarios. First, if the best motion 
vector candidate from step 110 is at the center of Vj (the "integer-pel vector") 130, then the 
encoder checks three candidates of 1/3-pel accuracy between the center vector and the 1/2- 

1 5 pel location with the next lowest RD cost 132. Second, if the best motion vector candidate 
from step 1 10 is a corner vector 134, then, the encoder checks the four vector candidates of 
1/3-pel accuracy that are closest to such corner 136. Third, if the best motion vector 
candidate from step 110 is between two corners 138, then, the encoder determines which of 
these two corners has lower RD cost and checks the four vector candidates of 1/3-pel 

20 accuracy that are closest to the line between such corner and the best candidate from step 1 1 0 
140. It should be noted that in implementing this process step 138 may be unnecessary 
because if V2 is neither at the center or a corner vector, then it would necessarily be between 
two corners. If the encoder is set to find motion vectors with 1/3-pixel accuracy, FIG. 7 
could be modified to end rather than continuing with step 114. 

25 

Computation And Memory Savings 

Because step 108 checks only motion vector candidates of 1/2-pixel 
accuracy, the computation and memory requirements for the hardware or software 
implementation are significantly reduced. To be specific, in a smart implementation 
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embodiment of this fast-search the reference frame is interpolated by 2x2 in order to obtain 
the RD costs for the 1/2-pel vector candidates. A significant amount of fast (or cache) 
memory for a hardware or software encoder is saved as compared to Telenor's approach that 
needed to interpolate the reference frame by 3x3. In comparison to the Telenor encoder, this 
5 is a cache memory savings of 9/4 or a factor of 2.25. The few additional interpolations can 
be done later on a block-by-block basis. 

Additionally, since the interpolations in step 108 are used to direct the search 
towards the lower values of the RD cost function, a complex filter is not needed for these 
interpolations. Accordingly, computation power may be saved by using a simple bilinear 

10 filter for step 108. 

Also, other key coding decisions such as selecting the mode of a macroblock 
(e.g., 16x16, four-8x8, etc.) can be done using the 1/2-pel vectors because such decisions do 
not benefit significantly from using higher accuracies. Then, the encoder can use a more 
complex cubic filter to interpolate the required sub-pixel values for the few additional vector 

15 candidates to check in the remaining steps. Since the macroblock mode has already been 
chosen, these final interpolations only need to be done for the chosen mode. 

Use of multiple- filters obtained computation savings of over twenty 
percent in running time on a Sparc Ultra 10 Workstation in comparison to Telenor's 
approach, which uses a cubic interpolation all the time. Additionally, the fast-memory 

20 requirements were reduced by nearly half. Also, there was little or no loss in 

compression performance. Comparing one preferred embodiment of the fast-search, 
Benzler's technique requires about 70 interpolations per pixel in the Telenor encoder and 
the present invention requires only about 7 interpolations per pixel. 

25 Coding The Motion Vector And Accuracies With Bits 

Once the best motion vector and accuracy are determined, the encoder 
encodes both the motion vector and accuracy values with bits. One approach is to encode 
the motion vector with a given accuracy (e.g., half-pixel accuracy) and then add some extra 
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bits for refining the vector to the higher motion accuracy. This is the strategy suggested by 
B. Girod, but it is sub-optimal in a rate-distortion sense. 

In one preferred embodiment of the present invention, the accuracy of the 
motion vector for a macroblock is first encoded using a simple code such as the one given 
5 in Table 1 . Any other table with code lengths {1, 2, 2} could be used as well. The bit rate 
could be further reduced using a typical DPCM approach. 



Code 


Motion 
Accuracy 


1 


1/2-pel 


01 


1/3-pel 


11 


1/6-pel 



Table 1, VLC table to indicate the accuracy 
mode for a given macroblock. 
Next, the value of the vector/s in the respective accuracy space is encoded. These bits can 

1 0 be obtained from entries of a single VLC table such as the one used in the H26L codec . The 
key idea is that these bits are interpreted differently depending on the motion accuracy for 
the macroblock. For example, if the motion accuracy is 1/3 and the code bits for the X 
component of the difference motion vector are 00001 1, the X component of the vector is 
Vx= 2/3. If the accuracy is 1/2, such code corresponds to Vx= 1. 

1 5 Compared to the Benzler method for encoding the motion vectors with a 

variable length code ("VLC") table that could be used for encoding 1/2 and 1/4 pixel 
accurate vectors, the method of the present invention can be used for encoding vectors of 
any motion accuracy and the table can be interpreted differently at each frame and 
macroblock. Further, the general method of the present invention can be used for any 

20 motion accuracy, not necessarily those that are multiples of each other or those that are of 
the type 1/n (with n an integer). The number of increments in the given sub-pixel space is 
simply counted and the bits in the associated entry of the table is used as the code. 



1 Observe that this code is the fourth entry (code number 3) of H26L's VLC table in [6], 
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From the decoder's viewpoint, once the motion accuracy is decoded, the 
motion vector can also be easily decoded. After that, the associated block in the previous 
frame is reconstructed using a typical 4-tap cubic interpolator. There is a different 4-tap 
filter for each motion accuracy. 
5 The AMA does not increase decoding complexity, because the number of 

operations needed to reconstruct the predicted block are the same, regardless of the 
motion accuracy. 

Experimental Results 

10 FIGS. 8-18 show test results of the Telenor encoder codec with and without 

AMA in a variety of video sequences, resolutions, and frame rates, as described in Table 2. 
These figures show rate-distortion ("RD") plots for each case. The "Anchor" curve shows 
RD points from optimized H.263+ (FIGS. 8 and 9 only). The "Telenor 1/2+b" curve shows 
Telenor with 1/2-pel vectors and bilinear interpolation (the "classical case"). The "Telenor 

15 1/3" curve shows the current Telenor proposal (the "Telenor encoder"). The 
"Telenor+AMA+c H curve shows the Telenor encoder with the full-search strategy of the 
present invention. The " Telenor +FSAMA+c", as shown in FIGS. 15-17, shows the current 
Telenor encoder with the fast-search strategy. (Unless otherwise specified, the full-search 
version of AMA was the encoder strategy used in the experiments.) All of the test results 

20 were cross-checked at the encoder and decoder. These results show that with AMA the gains 
in peak signal-to-noise ratio ("PSNR") can be as high as 1 dB over H26L, and even higher 
over the classical case. 



Video sequence 


FIG.# 


Resolution 


Frame rate 


Container 


FIG. 8 


QCIF 


10 


News 


FIG. 9 


QCIF 


10 


Mobile 


FIG. 10 


QCIF 


10 
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FIG. 11 


SIF 


15 


Garden 


FIG. 12 


QCIF 


15 


Tempete 


FIG. 13 


SIF 


15 




FIG. 14 


QCIF 


15 


Paris Shaked 


FIG. 15 


QCIF 


10 



Table 2. Description of the Experiments 
The video sequences are commonly used by the video coding community, 
except for "Paris Shaked." The latter is a synthetic sequence obtained by shifting the well- 
known sequence "Paris" by a motion vector whose X and Y components take a random value 
5 within [-1,1]. This synthetic sequence simulates small movements caused by a hand-held 
camera in a typical video phone scene. 

Comparison Of Full-Search And Fast-Search AMA 

The experimental results shown in FIGS. 16 and 17 demonstrate that the 
10 encoder performance with fast-search ("Telenor FSAMA+c") and full-search ("Telenor 
AMA+c") strategies for AMA is practically the same. This is true because the fast-search 
strategies exploit the convexity of the RD cost curve in the sub-pixel velocity space. In other 
words, since the shape of the RD cost follows a smooth convex curve, its minimum should 
be easy to find with some smart fast-search schemes that descend down the curve. 

15 

Combining AMA And Multiple Reference Frames 

In the plot shown in FIG. 18, the curves labeled "lr" used only one reference 
frame for the motion compensation, so these curves are the same as those presented in FIG. 
10. The curves labeled "5r" used five reference frames. 
20 The experiments show that the gains with AMA add to those obtained using 

multiple reference frames. The gain from AMA in the one-reference case can be measured 
by comparing the green and pink curves, and the gain in the five-reference case can be 
measured between the blue and red curves. 



It should be noted that the present invention may be implemented at the 
frame level so that different frames could use different motion accuracies, but within a frame 
all motion vectors would use the same accuracy. Preferably in this embodiment the motion 
vector accuracy would then be signaled only once at the frame layer. Experiments have 
5 shown that using the best, fixed motion accuracy for the whole frame should also produce 
compression gains as those presented here for the macroblock-adaptive case. 

In another frame-based embodiment the encoder could do motion 
compensation on the entire frame with the different vector accuracies and then select the best 
accuracy according to the RD criteria. This approach is not suitable for pipeline, one-pass 

10 encoders, but it could be appropriate for software-based or more complex encoders. Still 
another fame-based embodiment the encoder could use previous statistics and/or formulas 
to predict what will be the best accuracy for a given frame (e.g., the formulas in set forth in 
the Ribas work or a variation thereof can be used). This approach would be well-suited for 
one-pass encoders, although the performance gains would depend on the precision of the 

1 5 formulas used for the prediction. 

The terms and expressions which have been employed in the foregoing 
specification are used therein as terms of description and not of limitation, and there is no 
intention, in the use of such terms and expressions, of excluding equivalents of the 
features shown and described or portions thereof, it being recognized that the scope of the 

20 invention is defined and limited only by the claims that follow. 
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WHAT IS CLAIMED IS: 

1 . A fast-search adaptive motion accuracy search method for estimating 
motion vectors in motion-compensated video coding by finding a best motion vector for a 
macroblock, said method comprising the steps of: 

(a) searching a first set of motion vector candidates in a grid of sub-pixel 
resolution of a predetermined square radius centered on V t to find a 
best motion vector V 2 ; 

(b) searching a second set of motion vector candidates in a grid of sub- 
pixel resolution of a predetermined square radius centered on V 2 to 
find a best motion vector V 3 ; and 

(c) searching a third set of motion vector candidates in a grid of sub-pixel 
resolution of a predetermined square radius centered on V 3 to find 
said best motion vector of said macroblock. 

2. The method of claim 1, said step of searching a first set of motion 
vector candidates in a grid of sub-pixel resolution of a predetermined square radius centered 
on V x to find a best motion vector V 2 further comprising the step of searching a first set of 
eight motion vector candidates in a grid of 1/2-pixel resolution of square radius 1 centered 
on Y x to find a best motion vector V 2 . 

3. The method of claim 1, said step of searching a second set of motion 
vector candidates in a grid of sub-pixel resolution of a predetermined square radius centered 
on V 2 to find a best motion vector V 3 further comprising the step of searching a second set 
of eight motion vector candidates in a grid of 1/6-pixel resolution of square radius 1 centered 
on V 2 to find a best motion vector V 3 . 

4. The method of claim 1 further comprising the steps of using V 2 as the 
motion vector for the block if V 2 has the smallest rate-distortion cost and skipping step (c) 
of claim 1 . 
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5. The method of claim 1, said step of searching a third set of motion 
vector candidates in a grid of sub-pixel resolution of a predetermined square radius centered 
on V 3 to find said best motion vector of said macroblock further comprising the step of 

5 searching a third set of eight motion vector candidates in a grid of 1/6-pixel resolution of 
square radius 1 centered on V 3 to find said best motion vector of said macroblock. 

6. The method of claim 1, said step of searching a third set of motion 
vector candidates in a grid of sub-pixel resolution of a predetermined square radius centered 

10 on V 3 to find said best motion vector of said macroblock further comprising the step of 
skipping motion vector candidates of said third set of motion vector candidates that have 
already been tested. 

7. The method of claim 1 further wherein said step of searching said first 
15 set of motion vector candidates further comprises the step of searching said first set of 

motion vector candidates using a first filter to do a first interpolation, said step of searching 
said second set of motion vector candidates further comprises the step of searching said 
second set of motion vector candidates using a second filter to do a second interpolation, and 
said step of searching said third set of motion vector candidates further comprises the step 
20 of searching said third set of motion vector candidates using a third filter to do a third 
interpolation. 

8. The method of claim 1, said step of searching a second set of motion 
vector candidates in a grid of sub-pixel resolution of a predetermined square radius centered 

25 on V 2 to find a best motion vector V 3 further comprising the steps of: 

(a) searching three candidates of 1/3-pel accuracy V 2 and a 1/2-pel 
location with the next lowest RD cost if V 2 is at the center; 

(b) searching four vector candidates of 1/3-pel accuracy that are closest 
to V 2 if V 2 is a corner vector; and 



18 



(c) determining which of two corners has lower RD cost and searching 
four vector candidates of 1/3-pel accuracy that are closest to a line 
between said corner with lower RD cost, if V 2 is between two corners 
vectors, 

5 

9. An adaptive motion accuracy search method for estimating motion 
vectors in motion-compensated video coding by finding a best motion vector for a 
macroblock, said method comprising the steps of: 

(a) searching a first set of motion vector candidates in a grid centered on 
10 Vj to find a best motion vector V 2 using a first filter to do a first 

interpolation; 

(b) searching a second set of motion vector candidates in a grid centered 
on V 2 to find a best motion vector V 3 using a second filter to do a 
second interpolation; and 

1 5 (c) searching a third set of motion vector candidates in a grid centered on 

V 3 to find said best motion vector of said macroblock using a third 
filter to do a third interpolation. 



10. The method of claim 9 wherein said step of searching using a first 
20 filter to do a first interpolation further comprises using a simple filter to do a coarse 
interpolation. 



11. The method of claim 9 wherein said step of searching using a first 
filter to do a first interpolation further comprises using a simple filter to do a coarse 
25 interpolation and said step of searching using a second filter to do a second interpolation 
further comprises using a complex filter to do a fine interpolation. 
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12. The method of claim 1 1 wherein said step of searching using a third 
filter to do a third interpolation further comprises using a complex filter to do a fine 
interpolation. 



5 13. The method of claim 9 wherein said step of searching using a first 

filter to do a first interpolation further comprises using a bilinear filter to interpolate the 
reference frame by 2x2. 



14. The method of claim 9 wherein said step of searching using a first 
10 filter to do a first interpolation further comprises using a bilinear filter to interpolate the 
reference frame by 2x2 and said step of searching using a second filter to do a second 
interpolation further comprises using a cubic filter to do a fine interpolation. 



15. The method of claim 14 wherein said step of searching using a third 
15 filter to do a third interpolation further comprises using a cubic filter to do a fine 
interpolation. 



16. An adaptive motion accuracy search method for estimating motion 
vectors in motion-compensated video coding by finding a best motion vector for a 
20 macroblock, said method comprising the steps of: 

(a) searching at a first motion accuracy for a first best motion vector of 
said macroblock; 

(b) encoding said first best motion vector and said first motion accuracy; 

(c) searching for at least one second best motion vector of said 
25 macroblock at an at least one second motion accuracy; 

(d) encoding said at least one second best motion vector and said at least 
one second motion accuracy; and 

(e) selecting the best motion vector of said first and at least one best 
motion vectors using rate-distortion criteria. 
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1 7 . The method of claim 1 6 wherein said step of selecting the best motion 
vector using rate-distortion criteria further comprises the step of said rate-distortion criteria 
adapting according to the different motion accuracies to determine both the best motion 

5 vectors and the best motion accuracies. 

1 8 . The method of claim 1 6, said step of searching for at least one second 
best motion vector at an at least one second motion accuracy further comprising the step of 
searching for at least one second best motion vector of said macroblock at an at least one 

10 second motion accuracy that is finer than said first motion accuracy. 



vector using rate-distortion criteria further comprises the step of using rate-distortion criteria 
of the type "distortion + L*Bits" to select the best motion vector. 



20. An adaptive motion accuracy search method for estimating motion 
vectors in motion-compensated video coding by finding a best motion vector for a 
macroblock, said method comprising the steps of: 

(a) searching at a motion accuracy for a best motion vector of said 



19. The method of claim 1 6 wherein said step of selecting the best motion 
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macroblock; 



(b) 



encoding said motion accuracy using a code from a VLC table that is 
interpreted differently at different coding units according to the 
associated motion vector accuracy; and 



(c) 



encoding said best motion vector in the respective accuracy space. 
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21. A system for estimating motion vectors in motion-compensated video 
coding by finding a best motion vector for a macroblock, said system comprising: 
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(a) a first encoder for searching a first set of motion vector candidates in 
a grid of sub-pixel resolution of a predetermined square radius 
centered on V t to find a best motion vector V 2 ; 

(b) a second encoder for searching a second set of motion vector 
candidates in a grid of sub-pixel resolution of a predetermined square 
radius centered on V 2 to find a best motion vector V 3 ; and 

(c) a third encoder for searching a third set of motion vector candidates 
in a grid of sub-pixel resolution of a predetermined square radius 
centered on V 3 to find said best motion vector of said macroblock. 

22, The system of claim 21 wherein said first, second, and third encoders 
are a single encoder. 
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METHODS FOR MOTION ESTIMATION WITH 
ADAPTIVE MOTION ACCURACY 



ABSTRACT OF THE DISCLOSURE 

5 

Methods for motion estimation with adaptive motion accuracy of the 
present invention include several techniques for computing motion vectors of high pixel 
accuracy with a minor increase in computation. One technique uses fast-search strategies 
in sub-pixel space that smartly searches for the best motion vectors. An alternate 

10 technique estimates high-accurate motion vectors using different interpolation filters at 
different stages in order to reduce computational complexity. Yet another technique uses 
rate-distortion criteria that adapts according to the different motion accuracies to 
determine both the best motion vectors and the best motion accuracies. Still another 
technique uses a VLC table that is interpreted differently at different coding units, 

1 5 according to the associated motion vector accuracy. 
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FIGURE 1 



FIGURE 2 



The encoder searches 
for the best integer-pel vector V1 . 
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The encoder searches for 
the best 1/3-pixel accurate 
vector V1/3 near V1. 
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FIGURE 4 



The encoder searches 
for the best integer-pel vector V1 
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The encoder searches for 
the best 1/6-pixel accurate 
vector V 1/6 near V1. 



106 



o 


o 


o 


o o 


o 


o 


o 


o 


o 


o 


o 




o 








w 


o 


w 


P) 




n 


((Ilk 

Mil/ 




VIP 


p 


p 


V|P 


Pi 


p 

u 




w 








p 


V 


u 


V 


V 


o 


O 


o 


o o 


O 


o 


0 


O 


o 


o 


o 


• 


• 


• o 




o 


• 


• 


• 


o 


o 


o 


o 


O 0 


o 


o 


0 


o 


o 


o 


o 


• 


o 


• o 


• 


o 


• 


o 


• 


o 


o 


0 


• 


o o 


• 


o 


0 


dD 


0 


0 


o 


• 


0 


• 0 


• 


o 


• 


o 


• 


o 


o 


o 


o 


o o 


o 


o 


0 


o 


o 


o 



m 


1-pel 


• 


1/2-pel 


• 


1/3-pel 


o 


1/6-pel 



1/6 pixel 



V 



x 



FIGURE 3 
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FIGURE 6 



Checking the eight motion vector candidates in a grid of 
1/2-pixel resolution of radius 1 which is centered on V1. 
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Denoting the candidate that has the smallest RD cost as V2. 
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Checking the eight motion vector locations in a grid of 
1/6-pixel resolution of radius 1 that is now centered on V2. 



Yes 




h 
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Denoting the best motion vector of the eight as V3. 



h 

116 



Checking the new motion vector candidates in 
the grid of 1/6-pixel resolution of radius 1 that is 
centered on V3. Observe that some of the candidates in 
the grid have already been tested and can be skipped. 



h 
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Selecting the candidate with the smallest 
RD cost as the motion vector for the block. 



h 
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FIGURE 7 



110 from 
FIG. 7 
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Is V2 at the center? 



Yes 
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No 



Checking three candidates 
of 1/3-pel accuracy between the 
center vector and the 1/2-pel 
location with the next lowest RD cost. 



[s V2 a corner vector?. 



Yes 



No 
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Checking the four vector 
candidates of 1/3-pel accuracy 
that are closest to such corner. 




136 



Determining which of the two corners 

has lower RD cost and checking 
the four vector candidates of 1/3-pel 
accuracy that are closest to the line 
between such corner and V2. 
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As a named inventor, I hereby appoint the following registered practi tioners) to prosecute this app lication and to tra nsact ail business in the Paten t 
and Trademark Office connected therewith: Q customer Number [ 

OR 



Registered practitioners) name/registration number listed below 



Place Customer 
Number Bar Code 
label hnm 



Name 



Registration 
Number 



Name 



Registration 
Number 



Karen Dana Oster 
Erich W. Merrill, Jr. 
Bruce A, Kaser 



37,621 
31,982 
31,531 



Additional registered practitioners) named on supplemental Registered Practitioner Information sheet PTO/SB/02C attached hereto. 



Direct all correspondence to: □ Customer Number 

or Bar Code Label 



OR OS Correspondence address below 



Name 



Miller Nash LLP 



Address 



3500 U.S. Bancorp Tower 



Address 



111 S.W. Fifth Avenue 



City 



Portland 



State 



OR 



ZIP 



97204-3699 



Country 



USA 



Telephone 



(503) 224-5858 



Fax 



(503) 224-0155 



I hereby declare that all statements made herein of my own knowledge are true and that all statements made on information and belief are 
believed to be true; and further that these statements were made with the knowledge that willful false statements and the like so made are 
punishable by fine or imprisonment, or both, under 18 U.S.C. 1001 and that such willful false statements may jeopardize the validity of the 
application or any patent issued thereon. 



Name of Sole or First Inventor: 



d A petition has been filed for this unsigned inventor 



Given Name (first and middle tif any!) 



Family Name or Surname 



Jordi 



Ribas-Corbera 



Inventor's 
Signature 



Date 



Residence: City 



Redmond 



State 



WA 



Country 



USA 



Citizenship 



Spain 



Post Office Address 



6125 145th Court NE 



Post Office Address 



City 



Redmond 



State 



WA 



ZIP 



98052 



Country 



USA 



53 Additional inventors are being named on the _j__supplemental Additional lnventor(s) sheet(s) PTQ/SB/02A attached hereto 



[Page 2 of 2] 



PTO/SB/02A (3-97) 
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Under the Paperwork Reduction Act of 1995, no persons are required to respond to a collection of information unless it contains a 
valid OMB control number. 



Please type a plus sign (+) inside this box -> [ -f- j 



DECLARATION 



ADDITIONAL INVENTOR(S) 
Suppfemental Sheet 

Page J__ of _J_ 



Name of Additional Joint Inventor, if any: 



PI A petition has been filed for this unsigned inventor 



Given Name (first and middle [if any]) 



Family Name or Surname 



Jiandong 



Shen 



Inventor's 
Signature 




Date 




Residence; City 


San Jose 


State 


CA 


Country 


USA 


Citizenship 


China 


Post Office Address 


170 West Tasman 


Post Office Address 




City 


San Jose 


State 


CA 


ZIP 


95134 Countty USA 



Name of Additional Joint inventor, if any: 



| | A petition has been filed for this unsigned inventor 



Given Name (first and middie pf any]) 



Family Name or Surname 



Inventor's 
Signature 



Date 



Residence: City 



State 



Country 



Citizenship 



Post Office Address 



Post Office Address 



City 



State 



ZIP 



Country 



Name of Additional Joint Inventor, if any: 



f~| A petition has been filed for this unsigned inventor 



Given Name (first and middle [if any]) 



Family Name or Surname 



inventor's 
Signature 



Date 



Residence: City 



State 



Country 



Citizenship 



Post Office Address 



Post Office Address 



City 



jState 



ZIP 



Country 



+ 



Burden Hour Statement This form is estimated to take 0.4 hours to complete. Time will vary depending upon the needs of the individual case. Any 
comments on the amount of time you are required to complete this form should be sent to the Chief Information Officer, Patent and Trademark 
Office, Washington, DC 20231. DO NOT SEND FEES OR COMPLETED FORMS TO THIS ADDRESS. SEND TO: Assistant Commissioner for 
Patents, Washington, DC 20231 . 



